ON THE LIMIT OF LARGE SURFACE TENSION FOR A FLUID MOTION 
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. Abstract. We study the free boundary Euler equations in two spatial dimensions. We prove that 

if the boundary is sufficiently regular, then solutions of the free boundary fluid motion converge to 
solutions of the Euler equations in a fixed domain when the coefficient of surface tension tends to 
infinity. 
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1. Introduction. 

Consider the initial value problem for the motion of an incompressible inviscid fluid with free 
boundary, whose equations of motions are given by^ 

fj = —Vporj in (1.1a) 

div(n) = in r/(J7), (1.1b) 

P\d-n(v.) = ^■^ on dr]{^), (1.1c) 

r?(0) = id, ?7(0) = mq- (l-ld) 

where 17 is a domain in R"; r/(t, •) is, for each t, a volume preserving embedding r]{t) : — )• M"' 
representing the fluid motion, with t thought of as the time variable (physically r/(t, x) is the position 
at time t of the fluid particle that at time zero was at x); ' denotes derivative with respect to t; 
u : Q{t) — )• M" is a divergence free vector fleld on J7(t) defined by u = ?) o ri~^ and representing the 
fluid velocity, where ^}{t) = r](t){Q); A is the mean curvature of the boundary of the domain 0(t); p 
is a real valued function on Q{t) called the pressure; finally, is a non-negative constant known as 
the coefficient of surface tension, id denotes the identity map, uq is a given divergence free vector 
field on ri, and div means divergence. The unknowns are the fluid motion r] and the pressure p, but 



^Here x is the Lagrangian coordinate of the fluid particle whose Euler coordinate at time t is ri{t, x) 
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notice that the system (1.1) is coupled in a non-trivial fashion, in the sense that the other quantities 
appearing in (1.1), namely u, A and i^{t), depend explicitly or implicitly on ij and p. 

We shall prove that in two (spatial) dimensions, if the embeddings rj are sufficiently regular at the 
boundary of Q, the boundary at time zero has constant mean curvature, and uq is tangent to the 
boundary, then, as A; — oo, solutions {r],p) to (1.1) converge to solutions to the incompressible Euler 
equations on the fixed domain Q, given by 

'C = -VpoC, (1.2a) 

< div(CoC~^) = (1.2b) 

^ C(0) = id, C(0) = no (1.2c) 

where ({t, •) is, for each t, a volume preserving diffeomorphism ((t) : 17 — t- fi. It is well known that 
the pressure p in the incompressible Euler equations is not an independent quantity, being completely 
determined by the velocity vector field = C ° (see e.g. [22] or section 7). 

In order to explain the aforementioned regularity of rj at the boundary (whose seeming necessity 
is discussed further below) as well as to state the main result, we need to introduce some definitions. 

Given manifolds M and N, denote by H^{M, N) the space of maps of Sobolev class s between M 
and N. For s > § + 1 define 

£'^{n) = = jry G H'{Q,W) J{r]) = 1,7]'^ exists and belongs to if'(?7(17), 17)}, 

where J is the Jacobian. <f^(17) is therefore the space of if*-volume-preserving embeddings of 17 into 
W. Define also 

P^(17) =V'^ = ^r]e i?'(17,M") J{r]) = 1, ?? : 17 ^ 17 is bijective and t?"^ belongs to i?'}, 
so that P^(17) is the space of if'^-volume-preserving diffeomorphisms of 17. Notice that ©^(fi) C 

Let B^^"^ (dQ) be the open ball about zero of radius 6o inside i?*"'"^/^(517). We shall prove that if 
6q is sufficiently small, then the map 

m = f, 

where / satisfies 

J(id+V/) = 1 in 17, 
f = h on 917, 



is a well defined map, and y^{Bg^ ^ {d^)) is a smooth submanifold of H'^~^^{Q). From this we then 
construct another map 

: P^(17) X ^{Blf^dn)) ^ ^^(17), 
defined by $(/?, /) = (id +V/) o /3, 
and we define (f;(17) C .^^(17) by 

^;(17) = cl.(p^(17)x^(^;+^(ai7))). 

Notice that since /? € P^(17), we have /3(517) = 917. Therefore, solutions ?y to (1.1) which belong 
to (f^(17) decompose in a part fixing the boundary and a boundary oscillation. This decomposition 
is one of the main ingredients of our proof. We shall also show that, under our hypotheses, V/ is 
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in fact 11/2 degree smoother than rj (though /3 is as regular as rj). This is the sense in which we 
work with embeddings which have smoother boundary values. The reason why one needs to consider 
(5^(0) rather than £^{^) is explained in section 1.1. 
We are now ready to state our main result. 

Theorem 1.1. Let be a domain in with smooth boundary dfl, and assume that the mean 
curvature of dQ is constant. Given k G M+ and a divergence free vector field uq € H^^^^{Q,'M?) 
tangent to the boundary, let {r]k,Pk) G C^{[^,Tk),£fi ^^{^)) x ([0, T^), be a solution to 

(1.1) with initial condition uq, and C, € C"*^ ([0, cxo), "D^ ^^(fi)) be a solution to (1-2), also with initial 
condition uq; r]k, Pk end C exist and are unique by [7, 22], and C, is defined for all time by [18]. 

Assume further that rj^ S ^^(^)- 

Then there exist a T > and a /cq G M+, such that ifTk is maximal, then > T for all k > ko, 
and r]k{t) C(t) as a curve in S'^ ^^{^) as k ^ 00. 

The reason for requiring uq G H^^^^, whereas solutions r] are only in H^^^^, is to ensure uniqueness 
(see [7] ) . We assume that the initial velocity uq is the same for all k for simplicity. The result would 
still hold for appropriate sequences {u^k} which converge, in a precise sense, to uq when k ^ 00. 
Note that the assumption that dQ has constant curvature implies that it is connected, for if dO, 
had multiple components they would each have to be a circle and all would have the same radius. 
Since is a bounded domain in the plane, one of these boundary components would also bound 
the non-compact component of — $7 and this circle would have to have larger radius than the 
others. Hence this circle must be the only component, so di} must be connected. However, the 
assumption of constant mean curvature at time zero is necessary. We comment on these remarks, 
and other matters, in section 7.1. We also point out that, in general, the convergence %(t) — )• C(*) 
is not expected to hold in (see the analogous results in [16, 20]). 

Remark 1.2. Throughout the paper, we shall use the fact that the results in [7] also give rj G 

Notation 1.3. We reserve O for the fixed domain, with i}{t) being always the domain at time t, i.e., 
0(t) = ri{t){Q,). Of course, $1(0) = ft. In several parts of the paper the subscript k will be dropped 
for the sake of notational simplicity. 

The mathematical study of equations (1.1) has a long history, although for a long time only results 
under restrictive conditions had been achieved. In particular, a great deal of work has been devoted 
to irrotational flows, in which case the free boundary Euler equations reduce to the well-known 
water-waves equations. See [28, 33, 5, 1, 2, 32, 24]. 

Not surprisingly, when equations (1.1) are considered in full generality, well-posedness becomes a 
yet more delicate issue, and most of the results are quite recent. In this regard, Ebin has showed that 
the problem is ill-posed if A; = [14], although Lindblad proved well-posedness for A: = when the 
so-called "Taylor sign condition" holds [25] (the linearized problem was also investigated by Lindblad 
in [26]). When k > 0, a priori estimates have been obtained by Christodoulou and Lindblad [6] and 
Shatah and Zeng [31], with well-posedness being finally established by Coutand and Shkoller [7, 8] 
(see also [30]). Other recent results, including the study of the compressible free boundary Euler 
equations, are [9, 10, 12, 11]. 

1.1. The physical significance of theorem 1.1 and the boundary regularity of (1.1). The- 
orem 1.1 not only gives a satisfactory answer to the natural question of the dependence of solutions 
on the parameter k; it also addresses a well motivated problem in Applied Science, namely, when 
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one can, by considering a sufficiently higli surface tension, neglect the motion of the boundary in 
favor of the simpler description in terms of the equations within a fixed domain. 

The physical intuition behind theorem 1.1 is very simple, as we now explain. The system (3.2) 
can be derived from an action principle with Lagrangian 

Cirj)=K{rj)-V{rj), (1.5) 

where 

K{v,v) = l I \n? (1-6) 

is the kinetic energy and 

y(r/) = k\dVt{t)\ - k\dVt{^)\ = k(^Area{dn{t)) - Area(5Jl(0))) (1.7) 
is the potential energy^. 

Our theorem 1.1 is almost an example of a general theorem on motion with a strong constraining 
force [16]. For the general theorem we are given a Riemannian manifold M and a submanifold A^. 
Also given is a function V : M —^M which has A as a strict local minimum in the sense that W = 
on A and D^V is a positive definite bilinear form on the normal bundle of A in M. Then if rik{t) is a 
motion given by the Lagrangian C{r], rf) = |(?), rf) — kV{ri) where ( , ) is the Riemannian metric, and 
if C,{t) is a Lagrangian motion in A of ^{C;C) with the same initial conditions as r]k{t), the theorem 
says that r/fc(t) converges to C{t) as /c — )• oo. Also % — > (, but the second derivative in general does 

not converge. For our theorem, M = £^'^\n), A = vf,'^\n), ( , ) is the inner product on tangent 
vectors and V{r]) is given by (1.7). 

Our theorem 1.1 is not actually an example of the general theorem for two reasons: 

a) The inner product on tangent spaces is only a weak Riemannian metric. The topology that 

it induces is weaker than the H^^^^ topology of ^^{^)- 

b) The bi-linear form D^V is only weakly positive definite on each normal space; it gives the H^^^"^ 
topology rather than the H'^^^^ topology. 

We now turn our attention to the "boundary regularity" previously mentioned, i.e., to the seeming 
necessity of working with S'^{^) instead of As already pointed out, theorem 1.1 can be viewed 

as a generalization to infinite dimensions of established results about the behavior of the Euler- 
Lagrange equations near a submanifold which minimizes the potential energy. As in [16, 19, 20], in 
the present work the manifold minimizing the potential energy is 2?^ ''^(^)- Thus it is natural to try 
a similar approach to that of [16, 19, 20]. Let us outline what such an approach looks like, and then 
we shall point out why it is problematic, and how the introduction of ^^{^) solves the problem. We 
stress that the arguments sketched below are heuristic, their purpose being mainly a motivation for 
the definition of (^^'(J^). 

We would like to decompose the fluid motion as a part fixing the boundary and a boundary 
oscillation, and then produce estimates showing that the boundary oscillation goes to zero when 
k ^ oo. To this end, one first needs to show that is an infinite-dimensional Banach manifold 

with T>^{il.) as an embedded submanifold. Following well known arguments [15, 22, 13], £^{i^) can 



Many authors consider instead V{ri) = k\dQ,{t)\. As the equations of motion remain unchanged by adding a 
constant, here we choose to normalize the potential energy to make V = a.t time zero. Such a normalization is 
convenient for our purposes as we are interested in taking k oo, in which case, if we did not subtract the contribution 
at time zero, V{ri) would diverge to infinity. 
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then be endowed with a Riemannian metric which is inherited by D^(Q). As usual, one expects the 
resulting metric to be only a week Riemannian metric^. 
We then seek a decomposition of the form 

Ti,{£p = TpV^ + upiV^) = div-i(O), o p + o /3 (1.8) 

where /3 G 2?^ C f*, div~^(0)jy denotes divergence free vector fields which are tangent to the 
boundary, and I'lsiV^) denotes the normal space, at /3, of 'C^(i^) inside The next step would 

be to combine (1-8) with an implicit function type of argument in order to obtain 

f^«z.(P^) nearP^, (1.9) 

where i^i'Dfi) is the normal bundle of T>^{il.) considered as a submanifold of £"^(0). From these 
considerations it should follow that elements r] G £^{^) which are sufficiently close to P^(f]) can be 
written as 

r? = (id+V/)o/3, (1.10) 

giving the decomposition. 

Let us now see how the previous argument fails unless embeddings r]{il.) such that drji^D.) is suffi- 
ciently regular are considered. Suppose we succeed in constructing all the necessary manifold struc- 
tures, including the decomposition (1.8). To obtain the identification (1.9), we use a tubular neigh- 
borhood argument (see e.g. [23]), which in turn relies on of the exponential map: if V/ G z^id(I^^) is 
small enough, then tVf is a geodesic through the origin and then 

expid(V/) = (id+tV/) \t=i = id+V/, 

and at any other point /3 G P^(J7), 

exp^(V/) = (id+V/)o/3. 

Therefore, for t] G £^{^) sufficiently close to X'^(Q), we would obtain decomposition (1.10). 

Unfortunately the previous argument does not work. In our infinite dimensional setting, the 
existence of geodesies is not immediately guaranteed, since the metric is only a week Riemannian 
metric. In order to prove the existence of geodesies, one has to derive a well-posedness result for the 
following second order ODE on £^{Q), which is obtained as critical point of the kinetic energy for a 
curve e : [0,T] ^ <S^(0), 

(0 = (vAq 1 (div (Vs,e-^e o ) o in 0, 

1^(0) = id, ^(0) = V/, 

where Aq ^ means the inverse of the Laplacian (on the domain 9{Cl)) with zero boundary condition. 
To prove existence, we attempt a Picard iteration argument, which requires the right hand side of 
the equation to be a (sufficiently) smooth functions of its arguments, and this turns out not to be 
the case: eoO"^ e H"; V q^q-J o 9"^ G H""^; then div(Vg„g_i^ o 9"^) is also in H""^ since 9 o 9~^ is 

divergence free; but because 9 G H^, we have that 9 {qq G H'^^^ , i.e., d9{fl) is only H'^~^ regular and 
therefore ^(div(Vg„g_i^o0-i)) G from which follows that V(Ao^(div(V^^g_i^o0-i))) g H'-\ 



We recall that a weak Riemannian metric is one which induces, on each tangent space, a weaker topology than the 
one given by the local charts. This is a feature exclusive to infinite dimensional manifolds; see [15] for details. 
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i.e., there is a loss of a derivative. Notice that this has nothing to do with the regularity of d^l itself, 
as we have assumed that the domain $7 has smooth boundary^ ^. 

Here we avoid these problems by working explicitly in (^^^(0,), which has the decomposition (1.10) 
built into it. 

Notation 1.4. When the manifolds are clear from the context, we write simply or H^[M) for 
H^[M,N). Denote by || • \\s and || • \\s^d the Sobolev norms on Vt and dVt, respectively. In particular 
the norms are || • |[o and || • ||o,9. By Hq{Q) we shall denote the space of maps on Q, which 
vanish on dQ. 

1.2. Auxiliary results. Here we recall some well known facts which will be used throughout the 
paper. For their proof, see e.g. [15, 29, 4]. 

Proposition 1.5. Let s > ^ + 1, g e f G H'{n). Then foge H%n) and 

II /og||,<C 11/11,(1+ II 5 lis), 

where C = C{n, s, Q,). 

We shall make use of the following well-known bilinear inequality (see [15]) 

II UV \\r< II "U llrll V \\s, (l-H) 

for s > ^, s > r > 0, where C = C{n, s, r, Q). Recall also that restriction to the boundary gives rise 
to a bounded linear map, 

II u \\s,a< C II u ||^_^i, s > 0, (1.12) 

with C = C(n,s,0). 

Next we recall the decomposition of a vector field into its gradient and divergence free part. Given 
a vector field uj on 0, define the operator Q : H'{VL,W) ViJ'*+i(0, M") by Q{oj) = Vg, where 
g solves 

Ag = div{uj) in il., 
^ = (w, z/) on 557. 

Since solutions to the Neumann problem are unique up to additive constants, Vg is uniquely de- 
termined by so Q is weh defined. Define P : H'{n,W) div-i(O)^, where diy-^O)^ denotes 
divergence free vector fields tangent to dQ, hy P = I — Q, where I is the identity map. It follows 
then that Q and P are orthogonal projections in or that the image of Q is L^-orthogonal to the 
image of P. 

Notation 1.6. The letter C will be used to denote several different constants. C will never depend 
on k, rj or p. We will use both D and V to denote the derivative. 



^We notice that this lost of a derivative is an extra difficulty of the free boundary problem, not present in the case 
of the Euler equations in a fixed domain. In this situation, solutions are geodesies in T>^{Q,). Since the boundary 
is fixed, dil is, say, smooth (or at least assumed to be as regular as the rest of the data of the problem), and then 
Aq i(div(V9„e_i0 o 0-1)) e H'+\ so that V{Ao\dw{S/g^g-i6 o S'^))) G H". 

^Of course, showing that 9 is in the correct space is only a necessary condition to perform the Picard iteration, and 
does not by itself show that the ODE has a solution. The full proof has been carried out in the case of the Euler equations 
in a fixed domain by Ebin and Marsden [22], and more recently for geodesies on the group of symplectomorphisms by 
the second author [17]. 
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2. The space 

Here we construct the space as outhned in the introduction. We assume that s > f + 1, 

as usual. We shall make constant use of (1.11) to estimate the products involved. 
Start noticing that the condition J(id+V/) = 1 is equivalent to 

A/ + 0({D^ff + ■■■ + (D^fr) = 0. (2.1) 



Equation (2.1) can be considered as a non-linear Dirichlet problem for /, and so for / small, / should 
be determined by its boundary values. 

We shall present the argument for two dimensions, which is the main case of interest in this 
work. The interested reader can generalized the construction of <^^(f^) to higher dimensions. In two 
dimensions equation (2.1) becomes Af + fxxfyy — fxy = 0; or equivalently Af + det{D'^ f) = 0, where 
det means the determinant. Therefore, given we are interested in solving 

Af + det(L»V) = in 0, (2.2a) 
f = h on an. (2.2b) 



Define a map 



by F{h,f) = {f\on-h,Af + det{D^f)). 

Notice that F is in the neighborhood of the origin and and -F(0, 0) = 0, where we denote by 
the origin in the product Hilbert space H'^'^^{dQ) x H^~^{Q,). Letting w € H^~^^{Q,), we obtain 

D2FiO,0){w) = iw\9n,Aw), (2.3) 

where D2 is the partial derivative of F with respect to its second argument. From the uniqueness 
of solutions to the Dirichlet problem it follows that D2F{0, 0) is an isomorphism, and therefore by 

the implicit function theorem there exists a neighborhood of zero in //*^2 (90), which we can take 
without loss of generality as a ball '{dn), and a map ip : BI ^dn) H'+\n) satisfying 

V?(0) = 0, and such that F{h, (p{h)) = for all h € B^^ ^ {Q). In other words, / = ip{h) solves (2.2). 

Furthermore, since ip is continuous, given e > we can choose 5o so small that || (p{h) |[s+i< e. 
But if / is a solution of (2.2) with || / ||s+i< e, and e is sufficiently small, then by elliptic theory the 
solution is unique and obeys the estimate 

II / ||.+i< C II h l^^Q, (2.4) 

where the constant C depends only on e, s, and 

Next, we show that the map ip is injective. Let / = (p{h) and u = (p{g). Then f — u satisfies 




- det(L>2/) + det(L>2-u) in 
- g on 90, 

But 

det(i:'^u) - det(L>^/) = UxxUyy - uly - fxxfyy + fly 

— {Uxx fxx)Uyy + {fyy Uyy^fxx ~l~ {fxy ^xy^^fxy ~\~ Uxy^i 

and therefore combining (2.5) with (2.6) and standard elliptic estimates produces 

\\ f -u \\s+i< ci^W f -u \\s+i ( II / ||s+i + II u \\s+i )+\\h- g Ils+i,a (2-7) 



(2.5) 



(2.6) 
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where C depends only on and s. 

By (2.4) we can estimate / and u in terms of h and 5, respectively. Moreover, since || h q< 5q 

and II g W^j^i q< (5o, we can choose 5o so small that 

C \\ f -U \\s+l ( II / ||s+l + II ||s+l ) <]^\\ f - U \\s+l, 

and therefore this term can be absorbed on the left hand side of (2.7), yielding 

II / - ||s+i< C \\h- g lls+1,9 • 

In particular, f = u when h = g, so the map if is injective. 

Lastly, we show that the derivative of (p (which we already know exists and is continuous) is 
injective as well. Prom this and the above it then follows that (/^(fi^g (917)) is a submanifold of 
H'^~^^{Q). We start computing the derivative of ip at zero in the direction z. Since (p{0) = 0, 

ifiO + z)- ip{0) = ^{z) = Dp{0){z) + 0(11 z ll^^i^g). (2.8) 
But if{z) = w, where w solves 

Au; + det(i:»2^) = in 
w = z on dil, 

from which we can write 

w = -AQ\det{D^w)) + n{z), 

where Aq ^ is the inverse of the Laplacian with zero boundary condition and 'H(-z) is the harmonic 
extension of z to the domain Q,. From standard properties of the Laplacian and (2.4) we obtain 

II A^\det{D^w)) \Ui< C II z Wl^.^g . 

From this estimate and the expansion (2.8) we conclude that 

Dp{0) = n. 

In particular, the derivative of (p is injective at the origin. Since is C^, we conclude that Dip is 
one-to-one, provided 60 is taken sufficiently small. 

Recall now the definition (1.4). Notice that ^> is well defined (if 60 is small) and its image belongs 
to £^{^) since J(/3) = 1 and, by construction, J(id-|-V/) = 1. 

We have therefore proven: 

Proposition 2.1. Let s > ^ + 1 and let B^^ ^ [dO) he the open hall of radius 5q in H^~^2 (dQ). Then 

if 60 is sufficiently small, there exists an emhedding ip : B^^^{dQ) — > H^~^^{Q), given explicitly hy 
ip{h) = f, where f solves (2.2). Moreover, the map <I> given hy (1-4) is well defined. 

Definition 2.2. Under the hypotheses of proposition 2.1, we define <^^(^^) ^ f^(r2) by 

^;(j7) = ci>(p^(f])x^(s;;^(ai7))). 



in il, 


(3.2a) 


in r]{n) , 


(3.2b) 


on dT]{Q:), 


(3.2c) 


in r]{Q), 


(3.2d) 


on dT]{i}), 


(3.2e) 




(3.2f) 
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3. A NEW SYSTEM OF EQUATIONS. 

In this section, we shall derive a different set of equations for the free boundary problem (1.1). 
In light of proposition 2.1 and the hypotheses of theorem 1.1 we can, from now on, assume that 
solutions r] to (1.1) can be written as 

r/ = (id+V/)o/3, (3.1) 
with /3 e 2?^(f^), V/ G H^{il.) and with / satisfying (2.2a). We also observe that 

m = id, V/(0) = 0. 

It is customary to write the pressure as a sum of an interior and a boundary term, namely, 
P = Po + kAn, so that the system (1.1) takes the form 

' i] = —Vp 7] = — (Vpo + kVAn) ° V 
Apo = - div(V„ii) 

Po\dr,{n) = 
AAh = 
•^H\dn{n) = ^ 
I r?(0) = id, 77(0) = no. 

The energy for the fluid motion (3.2) is given by the sum of the kinetic and potential energies 
(1.6) and (1.7), respectively, 

= lljv\'+k[\dm\-m 

with r] : [0,T) £^{^l). This energy is conserved, and therefore 

m = l\\^o \\l (3.4) 

where we have used that fj = uo r] and r]{0) = id. 
Differentiating (3.1) in time, 

^ = {Vf + vDVf + v)op, (3.5) 

where v is defined by 

/3 = n o /3. (3.6) 
Using 7/(0) = id and 77(0) = uq, from (3.5) we obtain that 

^^o = V/(0)+^;o. 

Since uq is divergence free and tangent to the boundary, one sees that V/(0) = 0, and so 

Vq = Puo = Uq. 

Differentiating (3.5) again and using (3.2a) gives the following equation for V/ 

V/ + 2D^Vf + V/ + {v + v- Vv)D Vf + v + v-Vv = -Vp o (id +V/), (3.7) 
where the operator D^^, acting on a vector w, is given in coordinates by 

{Dl,wy = v^v'djiw\ (3.8) 



(3.3) 
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or in invariant form by 

Define L on the space of maps from to by L = id+D^/, Li = PL and L2 = QL, where P and 
Q have been defined in section 1.2. Notice that Li is invertible on the image of P if / is small, since 
in this case it will be close to the identity. Applying Q o (id— LL~^P) = Q — L2L^^P to (3.7) we 
obtain 

' Vf + 2{Q-L2L^^P)D,Vf 

+ {Q-L2Ll^P)Dl^Vf = -{Q-L2L^^P){Vpo(id+Vf)) in 0, (3.9a) 
i) + P{Vyv) + 2L^^PDyVf + L-{^PDl^Vf = in f], (3.9b) 

t V/(0) = V/(0) = 0, v{Q) = Puo = UQ. (3.9c) 
The system (3.9) will be used to derive estimates for V/ and v. 

4. Energy estimates. 

The main goal of this section is to derive energy estimates which will be used to control the 
boundary oscillation V/. 

Notation 4.1. From the assumptions that C has constant mean curvature, it follows that 17 
is a disk of fixed radius. In order to simplify the notation, we shall assume once and for all that the 
radius is one, therefore dO, = S^. Clearly there is no loss of generality in doing so. 

Taking into account the previous remark on notation, denote by r and i' the tangent and outer 
unit normal to = 517, so that r points counterclockwise-wise. One easily computes then 

and 

DrV = T. 

We shall also need the following. In two dimensions 



Df] 

so it follows that 



dxV^dyT]'^ - dyv^dxv'^ V -^^^^ ^^^^ y V -^^'^^ ^^'^^ 

since the Jacobian determinant of rj is one. Therefore, any Sobolev norm of {Dr])~^ is bounded by 
that of Drj, 

II (Drj)-^ \\s< C II D'n \\s . (4.1) 

A similar statement holds for (D/3)~^. Moreover, changing variables, using /3(r2) = 0, and J(/3) = 1 
gives 



Jn J/3o/3(n) Jn 

||.< C II llo +C II Dp-^ \U-i< C II (3 llo +C II D(3-^ 
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Using again change of variables and /3(0) = Q and J(/3) = 1, 

II Dp~^ ||,„i< C II (Z)/3)-i |l,„i 

so that 

II Ils< C II (3 \\s . (4.2) 

Given an initial condition uo S H^^^^{n) for the free boundary problem, by the results of Coutand 
and Shkoller [7], there exist a > 0, a unique curve T]k{t) G H'^^^'^^^l), and a unique function 
Pk{t) £ H'^{Q{t)) satisfying the free boundary value problem (3.2); moreover rjk ^ H^^/^{Q). 

Consider decomposition (3.1). As mentioned in the introduction, we shall prove that V/ is in fact 
11/2 degree smoother than rj; see proposition 4.4. Therefore it makes sense to consider V/^ G if^(r2), 
41/2 < s < 6. From (3.9c) we know that V/fc and V/fc are small for small time. More precisely we 
shall assume the following 

Assumption 4.1. (Bootstrap assumption). There exists constant C > 0, independent of k, such 
that on some interval [0, T^), 

II V/fc ||,< < s < 6, (4.3a) 

k 2 

II V/fc ||,< < s < 41/2. (4.3b) 

k 2 



Assumption 4.1 will be assumed in this section. Its validity is shown in proposition 4.4, which is 
stated below and proven in section 5. 

We seek to show that estimates (4.3) hold on a time interval uniform in k] in other words, that 
Tk does not become arbitrarily small as k gets large. In order to do this we shall prove that (4.3) 
implies improved estimates for V/ and V/, possibly with a different constant C > 0. 

Remark 4.2. For simplicity, in (4.3a) and (4.3b), all values < s < 6 and < s < 41/2 are 
considered. As we point out in section 6, it will be enough to have (4.3) for only finitely many values 
of s. 

Remark 4.3. Several estimates below are valid for || V/fe ||s and || V/fc \\s sufficiently small. In light 
of assumption 4.1, and since we are interested in the limit k — )• 00, there will be no loss of generality in 
assuming that such smallness conditions are in fact met. The most common application of this idea 
will be to drop higher powers of || V/^ ||s and || V/^ ||s, and also to keep only the highest derivative 
terms. For example, if one has a long expression such as C \\ V fk Us (1+ || V/^ II4) + C || V/fc 1^ 
+C II V/fc II21/2' if ^ is large enough this expression is controlled by C(|| V/fc II3 + || V/fc ||2i/2)- 
When the norms dealt with are clear from the context, we shall use a slight abuse of language and 
refer simply to V/fc or V/fc as being small. 

The propositions below will be used in obtaining the desired energy estimates. Their proofs are 
given in section 5. 

Proposition 4.4. Let r/fc S ^'^{^) be a solution of the free boundary problem (3.2) as in theorem 
1.1, and defined on a time interval [0,Tfc). Then rjk G H'^^/'^{Q). Moreover, ifT^ is sufficiently small 
and k sufficiently large, then V/fc € H^^i}), V fk £ H^^/'^iQ), and equations (4-3) hold. 



Proof. See section 5. 



□ 
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It will also be important to relate equations in ^1 and Q(t). To this end, let ^ denote either r] or 
id +V/ = r] o I3~^. We denote by the operator defined by 

A^g = {A{g o C')) o 

where A is the Laplacian in the domain ^(^2). 

Proposition 4.5. For sufficiently small time and small V/, A^ is an invertible elliptic operator on 

Proof. See section 5. □ 

We are now ready to start our energy estimates. To facilitate the reading, most of the calculations 
will be split into several short lemmas, propositions and corollaries. 
Define 

^ = id +V/. 

By construction, rj is volume preserving, so the vector field 

u = rj o rj~^ (4.4) 

is divergence free. Right invariance through f3 implies that rj{^) = drj{0,) = dr]{^}). In 

particular, p is defined on rj{i^)- 
Define the energy 

m = l I \^^+k\d^{si)\ 

^ (4.5) 
+ k\dr]{n)\. 



n 



The purpose of using E rather than the physical energy defined in (3.3) is that E allows us to 
estimate / more explicitly. Notice, however, that as E does not carry direct physical significance, 
we have not normalized it as to keep it finite in the limit k ^ oo. 

Lemma 4.6. 

m = h\^o iio +km + f [ + f [ vpo?j). (4.6) 

^ Jo Jn Jo Jn 

Proof. Direct calculation gives, 

— \drj{n)\= / Ao?j{?j,N o?j)\Drrj\, 
J on 

where N is the normal to dTj(Q), Dr is the derivative tangential to dVt, and \Drrj\ is the length of 
the derivative of Tj applied to the boundary surface element. Consider 

{Vp,rj orj~^) = — / pdiv(?7 o ry""*^) + / p{rjorj~^,N) 
fj{n) Jrj{n) Jdriin) 

p o rj {rj, N Off) \D.,-r]\ = k / A o rj {rj, N orf) \D.,-r]\, 
dn Jdn 

where we have used div {rj o rj~^) = and porj= kA orj on d^. Hence 

d_ 



k-^\dr]{Q)\ = I {Vp,i]or] ^) = / (Vp 077,77). 
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Now the result is straightforward from the fundamental theorem of calculus and the definition of 
E. □ 

Lemma 4.7. 

Jn Jn Jn 

Proof. Since rj = V/, recalling that L2 = QL, L = id +D^f and using the fact that Q is symmetric 
with respect to the inner product, 

/ {Q-L2L-'P){Vpo^)) = [ (V/, (id-(id+Z)2/)^r'^)(Vpo^)) 
Jn Jn 

= [ (V/, VpoJf)- [ (V/, D'fL^'P{Vpo7j)), 
Jn Jn 

where we used that Q and P are orthogonal and L^^ takes the image of P onto itself. □ 
Lemma 4.8. 



E{t) = k\dn\+[ [ {Vf,D^fL^'P{\/po?j))- [ [{Vf,DlVf) 
Jo Jn Jo Jn 

+ 2 / / (V/, D^fL^^PD^Vf) + [ [ (V/, D^fL^^PDl.Vf). 
Jo Jn Jo Jn 



'0 Jn Jo Jn 

Proof. Notice that uq = because of (3.9c) and (4.4). Recall that rj = V/ and apply lemmas 4.6 
and 4.7 to get 

ft r ft 



Eit) = k\dn\+ [ [ {Vf,D^fL-^P{Vpo?j))+ [ /(V/, V/ + (Q-L2L-ip)(Vpo^)). 
Jo Jn Jo Jn 

We claim that 

/ (V/, D,Vf) = 0. (4.7) 
Jn 

In fact, integrating by parts and using div(t;) = produces 

/ (V/, D,Vf) = - [ {D,Vf, V/) + / (V/, Vf){v,v) = - [ {D„Vf, V/), 
Jn Jn Jdn Jn 

since (i^, f ) = 0. Using the form of L2 and L, from (4.7), equation (3.9a), and the properties of Q 
and P (as in lemma 4.7), we get the result. □ 

Lemma 4.9. 

^ II V/ ||2 +\dn\ + \l [\DrDf\^ - (r, DrDff) < ^E{t) - j (r, D^Df) + I \D^Df\^ 

+ 1 [ \DTDf\''{T,DrDf)+ [ R, 

^ Jdn Jdn 



where 



R = ^(^\DrDf\^ + 2(r, DrDf)f j\l - tf (l + t {\D^Df\^ + 2(r, D^Df)) ) ' dt. 
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Proof. Making a change of variables in (4.5) one has 



\dn{t)\ = Area(r?(aO)) = Aiea(jj{dn)) = [ \DrT] 

Jan 



The energy can then be written 



But 



E{t) = ^ [ \^\^ + k [ \DrV\ = l\\Vf\\l+k[ \DrV\. (4.8) 
^ Jn Jan ^ Jan 

[ \D^v\ = I \T + DrDf\, (4.9) 
Jan Jan 



and direct computation gives 



|r + DrDfl = ^/{T + DrDf,T + DrDf) 



= V^ + \DrDf\^ + 2{r,DrDf), 
where we used (t, r) = 1. Using Taylor's theorem with integral remainder gives 

y^l + \DrDf\^ + 2{T,DrDf) = 1 + ^{\DrDf\^ + 2(t, DrDf) 

-l(\DrDff + 2{T,DrDf)]\R, 



(4.10) 



with R as in the statement of the lemma. From this, (4.8), (4.9) and (4.10) the result follows. □ 
Lemma 4.10. 

/ (DrD.ff <C [ (Drff + C [ {{DrDfl'' - (t, D^Df)') (4.11) 
Jan Jan Jan 

Proof. Writing 

Df = DrfT + D,fu (4.12) 

we obtain 

DrDf = [Dlf + D,f)T + {DrD.f - Drf)v, (4.13) 

from which it follows that 

\DrDf\' = {Dlf + D,ff + [DrD.f - Drff, (4.14) 

and 

{T,DrDf) = Dlf + D,f. (4.15) 

Using (4.14) and (4.15) one obtains 

/ {\DrDf\'' - (r, DrDff) = [ {DrD.f - Drff. 

Jan Jan 

2 1 2 

The Cauchy inequality with e, i.e., ab < ^ + e > 0, gives 

/ {DrD.f - Drff = [ {{DrD.ff + {Drff " 2DrDJDrf) 

Jan Jan 

> [ {{DrD.ff + {Drff -e{DrD,ff- -{Drff), 

Jan e 
and the result follows by choosing e sufficiently small. □ 
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Lemma 4.11. 



^ II V/ \\l +m + C II DrD.f \\lg< ^E{t) 

+C (ll Drf + II Df \\l + II Df \\, g (1+ II Df lla^^^g) II D^Df 



Proof. Given any e > 0, by the Sobolev embedding theorem, there exists a constant C > 0, depending 
only on dO, and e, such that 



DrDf |bo(oc)< C II DrDf h_^^^g< C II Df \\s 



Use (4.16) to find 



/ iDrDfl" < sup \DrDf\^ [ \DrDf\^ <C\\Df \\l \\ D^Df W^q 



(4.16) 



(4.17) 



Analogously, 



\D,Df\'{T,D,Df) < \DrDf\' <sMDrDf\ \D^DfY 
an Jan an Jan 



(4.18) 



<C\\Df Ws^^J D^Df 
Next, we use integration by parts to estimate the term 

/ {T,D,Df) = - [ {DrT,Df). 

Jan Jan 

Since the boundary is S^, we see that DrT = —i^, u being the outer unit normal, and hence the 
integral becomes 

/ {T,DrDf)= I {u,Df)= [ DJ. 

Jan Jan Jan 

Now we will use the equation for /. Integrate (2.2a) over Q, and integrate the Laplacian term by 
parts to find 

/ D,f=[ {fly - f^^fyy) <C\\D^f\\l<C II Df \\l . 

Jan Jn 



Therefore 



/ {r,DrDf) 
Jan 



<C\\Df 



(4.19) 



Now we estimate the remainder term from lemma 4.9. By inequality (4.16) and the bootstrap 
assumptions we obtain that 



sup sup 



te[o,i] an (i^t{\DrDf\'^ + 2{T,DrDf)) 



(4.20) 



16 DISCONZI AND EBIN 

provided that k is sufficiently large. Using {t^Dt-DJ) < |Z)^L)/| (since |r| = 1), for any point on $7 
we find that 

36 4 
(|Z).Z)/|2 + 2(r, D^Df)) <CY, \DrDff = CY^ \D^Df\'+' 

e=3 £=1 

< C^\D^Df\^ sup \D^Dff 

Hence from inequalities (4.20) and (4.21), and invoking (4.16) once more, we find that the remainder 
term R in lemma 4.9 obeys obeys the estimate 



/ R<C\\Df Is _,^^g\\D,Df Wis, 
Jan 2 



(4.22) 



provided that V/ is sufficiently small, i.e., that k is large enough. Using lemmas 4.9 and 4.10 along 
with (4.17), (4.18) and (4.22), produces the result. □ 

Lemma 4.12. L~^P is a bounded linear operator in , s > 0, ifVf is sufficiently small. 

Proof. Since L = id +D^f, for V/ small this operator is bounded, with bounded inverse on the 
image of P. P is an orthogonal projection in L^(fi) and it is also bounded in H^{Q), s > 0. □ 

Lemma 4.13. 

^ II V/ \\l + II DrD^f \\Iq< c( II Drf WIq + II Df \\l (4.23) 
+ 11 Df \\s^^^, (1+ II Df lla^.,,^) II DrDf Wig 

II V/||o||I?.'.V/||o+^^ II V/||o||D/||2+.||ir'^ll (11 Vpo^llo 

+ II D,Vf llo + II DlVf llo 

where || L^^P || is the operator norm of L^^P. 
Proof. From lemmas 4.8 and 4.11 it follows that 

^ II V/ Wl + II DrD^f Wlo < C[ II Drf Wld + II Df Wl 

+ II Df II 

(^+ II ll§+£,a) II DrDf Wo,d 



+ (W, D^fL^'PVpo^- ^ r f {Vf,Dl,Vf) 

+ (V/, D^fL^^PD^Vf) + jf[ (V/, D^fL^'PDlVf). 



Now the lemma follows from the Cauchy-Schwarz inequality, (1.11) and the fact that Q, and hence 
P and are bounded operators in L? (see lemma 4.12). □ 

Lemma AAA. 

II Df lU <C(|| Df W2+eW Df lU + II DJ Wi,d 
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Proof. Since our goal is to estimate V/, it can be assumed that Jqq / = 0. Then one has the following 
elliptic estimate 

II / II, < C( II A/ ||,_2 + II D,f IL_|,a )• (4.24) 
A/ 111 < C II {D^ff lU < C II D^f |bo(n)|| |U < C \\ D^f ||i+,|| D^' ||i 

2 2 ^ ' 2 2 



Use (2.2a) and the Sobolev embedding theorem to obtain 
< C II {D^ff 111 < ( 
<C\\Df h+eW Df lU 



so that 

ll/ll|<C'(|| Df h+eW Df III + 11 D,f ||i,a^ 
from which the result follows. □ 

Proposition 4.15. Let A be the mean curvature of the boundary of the domain r]{Q). Then ifS/f 
is sufficiently small, 

M o ^ - 1 ||^_^ , < C (ll DrDf W,^,^, + II DlDf W,^.^, 
where s' = max(s,e), and e is a fixed small number. 

Proof. Notice that the left hand side of the above inequality is well defined since r]{Vt) = rj{^)- The 
mean curvature at rj{^) is given by 

\Drri\ \\Drr]\ 

^ Dl?, + -^AD%Dr^Dr?j, 



|D^^|2 ^' \D^^i 
where is the normal to drj(Q). Since rj = id +Df, one has 

= T + D^Df, 

and 

Dl?}=-iy + DlDf. 

Notice that by (4.16), {Drrjl^^ is well defined for large k, and the same holds for the several expressions 
below derived from \DT-rj\~^. Then 

{D%DrV) = {-y,DrDf) + {T,DlDf) + {DlDf,DrDf), 



so that 



^°^N°^=J^{^-DlDf) 

+ ((-'^' DrDf) + (r, DlDf) + {DlDf, DrDf) ) (r + DrDf) 



(4.25) 



We compute 

\Dr?}\^ = {t + DrDf, T + DrDf) = 1 + 2{DrDf, t) + \D^Df\^ = 1 + Mo, 

where 

Mo = 2{DrDf, t) + \DrDf\\ (4.26) 
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But, 



\DrW 1 + Mo 



1 + Mi, 



where 



Since 



Ml = -Mo + Mg 



1-t 
(1 + tMoY 



dt. 



(4.27) 
(4.28) 



= {l + MiY, we have 
A oT] N o?j = V + Miu - {I + Mi)dId f 

+ (1 + Mi)2 ((-z., DrDf) + (r, 1)2/)/) + {DlDf, DrDf)) (r + I),Z)/) 



where 

From this it follows 

where 

and 

We now have 
where 



= 1/ + Mil/ - (1 + Mi)DlDf + (1 + MifNhiT + D,D/), 
M2 = {-u,DrDf) + (r,D2^/) + {DlDf,DrDf). 
\Ao^'^ = l + 2(1/, M3) + M| = 1 + M4, 
M3 = Mil/ + (1 + Mi)DlDf + (1 + Mi)2M2(r + DrDf), 
M4 = 2(i/,M3) +M|. 



(4.29) 



(4.30) 
(4.31) 



1^ o ?7| = VI + M4 = 1 + M5, 



1 



M5 = -M4--M| 

2 4 7o (l + tM4)2 



1-t 



dt. 



(4.32) 
(4.33) 



Moreover, ^ o ry > for large fc, so we can drop the absolute value and write 

^o^ - 1 = M5. 

Our goal now is to use the above expression to estimate || Aorj —\ q- We shall use (1.11) to 
estimate the several products involved. If s > 0, then the condition s + i>|^ = iis satisfied; if 
s = then we use || • ||i q< C \\ ■ \\i_^_^q and apply (1.11) to the ^ + e norm. Henceforth we write 
s', with s' defined as in the statement of the proposition, covering both situations. 
From (4.26), 

|2 



II Mo |L,+ i_,< C [\\ DrDf + II DrDf ||f,^,_^j < C II DrDf „ (4.34) 

where the last inequality holds if V/ is sufficiently small. Next, from (4.28), 



Ml 11,,+ i^g < C II Mo |L,+ i,, + II Mo 11^,^,^^ 



1 - t 



(1 + iMo)' 



dt 



< C II NL 



\\s' + \,d 



1 + 



1 - 1 



(4.35) 



(l + iMo)2 



dt 



LARGE SURFACE TENSION 19 

where the last inequahty holds if V/ is sufficiently small because of (4.34). The Sobolev norm of the 
integrand can be estimated in terms of the Sobolev norm of Mq and the Sobolev norm of the real 
valued function ht{x) = jj^^^^i < t < 1, as follows. For 6 > very small it holds that 

\\ht\\s,[^s,s]<C, (4.36) 

where || ht \\s^[-5^s] is the Sobolev norm of ht over the interval [—(5,5] and the constant C does not 
depend on t. But by (4.34), we have that ht o Mq : 50 — > M satisfies ht o Mo{x) G {—S,S) for all 
X G dQ (provided that k is large enough), hence (4.35) gives 

II Ml 11^,^1^^ < C II Mo (l + C II Mo ) < C II Mo W.^^g 

< II DrDf 11,,+ .^,, 

where (4.34) has been employed. To estimate M2, use (4.29) to get 

II M2 ||,,+ i g< C (11 D^Df + II D^Df ||,,^i 3) , (4.38) 



where we ignored higher powers of the terms involved since we can assume that k is sufficiently large 
and hence V/ sufficiently small. 
Prom (4.30), 

II Ms 11,,+ .,^ < C II Ml 11,,+ .^, +C[l+ II Ah IL,+ i,a ) II D^Df l,^.^^g 

+ C(l+ II Ah 1,^,9 y\\M2 h^,9 (1+ II DrDf l,^.^g ) (4.39) 

2 i 



<C(|| D^Df |L,+ i,, + || D^^Df ||,,^i_. 



where we have used (4.38) and (4.37) and again dropped higher powers of the terms involved. From 
(4.31) and (4.39) it then follows that 



II M4 ||,,^_i,3< C [\\ DrDf l,^^^g + II D^Df l,^.^g ) . (4.40) 

Finally AI^ is estimated from (4.32) with the help of (4.40); the integral term in (4.32) is estimated 
by an argument similar to that used to obtain (4.37); it yields 

II Ms ||,,+ i,3< C (11 DrDf l,^^^g + II D^Df l,^.^g) . (4.41) 

Combining (4.33) and (4.41) the result follows. □ 

Proposition 4.16. Let Ah be the harmonic extension of the mean curvature to the domain 
Then if V/ is sufficiently small, 



II S/Ah o ^ II. < C II {DTjr^ \\r (II DrDf \\^,^^g + II D';Df \1,+ i^q) , (4.42) 
where r = max(s, 1 + e), s' = max(s, e), and e is a fixed small number. 

Proof. Notice that the left hand side of (4.42) is well defined since r]{Q) = ri{Q). Consider the 
function Ah : ^2 — >■ M defined by Ah = Ah orj — 1. It satisfies 

(AfjAH = inn, 
\ Ah = Aorj — 1 on dQ. 
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Intuitively, Ah is the harmonic extension of the mean curvature subtracted off the contribution from 
the domain at time zero, from the point of view of the fixed domain il.. Invoking proposition 4.5 and 
standard elhptic estimates, we have 



\\AH\\s<C\\Ao^-l\\^_i^g. (4.43) 

Computing, 

VAh = V {Ah o?j-l) = \/Ah o rj Drj, 

hence, 

VAh =VAh {Dlf)"^. 

Using (1.11), (4.43) impHes 

II VAno'nWs < C II VAh \\s\\ {D7j)-^ \\r 

<\\Ah \\s+i\\iDvr^ \\r (4.44) 
<|| ^o^-l ll^^i^gll (D^-i ||„ 

where r is as stated in the proposition. Combining (4.44) with proposition 4.15 gives the result. □ 
Proposition 4.17. 

II V \\s< C, 

for small time and s < 3. 
Proof. Define 

E, = ^\\v\\l=^{v,v)s, (4.45) 
where (•, ■)s is the inner product in H'^{0,). Differentiating in time and using (3.9b) produces 

= {v,v), = -{v,P{V,v))s - {v, 2L-^PV,Vf)s - {v,L^^PDl,Vf),. 
The last two terms obey 

{v, 2L-1PV,V/), <C\\v WlW V/ ||,+i 

(4.46) 

< CE^ II V/ ||s+i 

and 

{v,L^^PDl,Vf)s< C\\v\\l\\ V/ IU+2 

(4.47) 

< CE'S II V/ ||,+2 • 

In order to obtain a differential inequality for E^, we need now to show that (t',P(V„f))s can be 
estimated in terms of the s-norm of v. 

One can assume that s > 1 and also that s is an integer. In fact, when one or both of these 
conditions are not satisfied, we can replace s by the next larger integer and carry the estimates in 
this case, provided the corresponding norms of V/ and V/ are still bounded, which will be the case 
for s < 3 due to the bootstrap assumption. The condition s > 1 is required to use (1.11), whereas 
an integer s simplifies estimates with the Sobolev norms, as 



{v,P{V,v))s = Y,iD'v,D'P{V,v))o, 



£=0 
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where (•, •)o is the Lp' inner product. Since P is a bounded operator in the Sobolev norm, 
{D^v,D^P{V^v))q <|| D^v IIoII D^P{V^v) ||o< || v M P{V^v) \U 

< C II -u ll^ll V^v \\i< C II V \\}\\ V \\e+i< C II V |||+i, 

and so 

s 

I Y^iD'v, D'P{V,v))o < \{D'v, D'P{Vyv))o\ + C \\ v \\l . (4.48) 



1=0 



Consider 



{D'v,D'P{V,v))o= I {D'v,D'P{V,v)). 

Since P = I — Q we get 

/ {D'v,D'P{V^v)) = [ {D'v,D%V^v))+ [ {D'v, D'Q{V^v)). (4.49) 
Jn Jn Jn 



But 



[ {D'v,D'{V,v))= [ {D'v,V,{D%)) + [ {D'v,[D\VM 
Jn jQ Jn 

The first term on the right hand side is zero because integrating by parts and using {v, u) = gives 

{D'v,V,{D'v)) = - [ {V,{D'v),D'v). 



For the second term, notice that the commutator term involves at most s derivatives of v and hence 

3 



{D'v,D%Vyv)) = / {D'v,[D\V.,]v) <C\\v ||^< CE|. (4.50) 
! Jn 

Now consider J^{D^v, D'^Q('Vi,v)). By definition Q{Vvv) = V(7, where g solves 

{Ag = div(V„t;) in il, 
^ = {V,v,u) ondn. 

We observe the simple but important fact that since v is divergence free, 

div(Vt,t;) = diV^djV^, 

i.e., like in the ordinary case of the Euler equations, div(Vi,u) is only first order in derivatives of v 
rather than second order. Elliptic theory gives 

II Vfif \\s <\\ g \\s+i< C {\\ div^djv' \\s--i + || (V.„i;,z^) \\s-i,d) 

< C (ll Dv \\U + II (V,^;, u) ||^_ 1 < C (ll His + II (V.^;, z/) ||,_ i . 
To estimate the boundary term, recall that {v, v) = 0, so that (VuV, v) = —{v, Vt,J^) and hence 

II (V„«,z>) ||,_i,a=ll (^^, V„i/) \\,_i^g< C II V \\l_i^9< C II V \\l . (4.51) 
We conclude that || V^f ||s< C \\ v ||^, and therefore 

/ {D'v,D'Q{V^v)) = [ {D'v,D'Vg) <\\ D'v ||o|| I^'V^ 
Jn Jn 



<\\ V 11,11 II, < C II V f < CES. 



3 



(4.52) 
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Using (4.46), (4.47), (4.48), (4.49), (4.50) and (4.52) into (4.45) yields 

3 

Ev < C + C Ey , 



where we have used the hypotheses on V/ and V/. A standard ODE argument now shows that for 
small time E^ < C, what gives the desired result. □ 

Corollary 4.18. 

II /3 IU< C, 

for small time and s < 3. 

Proof. Since j3 = v o j3 and /3(0) = id, one has 

I3 = id+ f vop, 







and hence 

ft ft 



P II. < C + C I II V II, (1+ II p \\s) <C + Ct + C I 
Jo Jo 



where proposition 4.17, proposition 1.5 and J(/3) = 1 have been employed. Iterating the above 
inequality yields the result. □ 

Corollary 4.19. 

II ri \\s<C, 
II V \\s<C, 

for small time and s < 3. 

Proof. As before, we can assume that s > 2. Use then (3.5), proposition 1.5 and J{f3) = 1 to find 

II V \\s <|| V/ o /3 II, + II DyVf ol3\\, + \\vof3 II, 

< C II V/ II, (1+ II f3 II,) +C II V 11,11 D^fo/3\\s +C II V II, (1+ II /3 II,) 

<C{1+ II /3 II,) (ll V II, (1+ II V/ ||,+i)+ II V/ I 

Proposition 4.17, corollary 4.18 and the bootstrap assumption give || ?) ||,< C. The bound on ?y then 
follows from 7/ = id + Jq ?). □ 

Proposition 4.20. 

II Vpo o 1] \\s< C, 

for small time and s < 3. 

Proof. Define qo = Po o r] and notice that Vpo °f] = VqQ{Dri)^^. Therefore 

II Vpo ° V \\s< C II (Di])'^ y\\ Vqo ||„ (4.53) 

where s' = max{s, 1 + e}. We claim that both terms on the right hand side can be estimated from 
•fj. In light of (4.1), the first term is estimated in terms of || Drj ||,/, and hence, in terms of || r] ||,'+i. 
This last term, in turn, is estimated in terms of fj due to 

r? = id + / f], (4.54) 
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SO 



II {Dii)-^ h<C{l+ f W-nWs'+i). (4.55) 
For the second term on the right hand side of (4.53), notice that satisfies 

{^riflo = —{diU^djUi) o 1] in ft, 
qo = on 

so that 

II Vijo ||s<|| qo \\s+i < C II (diU^djUi) o T] \\s-i< C \\Vuor] \\1„ 
where we let s" = max{s — 1, 1 + e} and used (1.11). Since Vn o r] = V(n o r]){Dr])~^ , this imphes 
II Vgo II. < C II V{uo7]){D7])-^ \\l,< C II V(uor?) ||2„|| [Dr])-^ 

<C \\uor] ||^„+i (1 + /" Wv \\s"+i) 

Jo 

<C\\i) ||^„+i (1 + / II ?7 II," 
Jo 



2 

(4.56) 



^2 



after invoking (1.11), (4.1) and (4.54), and recalling that u o rj = rj. Using (4.55) and (4.56) with 
(4.53) yields 



II Vpo ° ^ lis < C II 7) ||2,^i (1 + / II 7] ||,' + l)(l + / II 7? ||s" + l)^ (4.57) 

Jo Jo 

so it is enough to estimate ?). 

Letting r = max{s" + 1, s' + 1, 2 + e}, we see that r < 3 and hence the result follows from corollary 
4.19. □ 

Proposition 4.21. 

II Vp o ^ ||o< C (1 + A: II V/ II2+. +k II V/ ||3+e) . 



Proof. From p = po + kAn one obtains 

II Vpo ^ ||o<|| Vpo o ^ llo +k II VAh °v\\o ■ 
Propositions 1.5 and 4.20, corollary 4.18 and J(/3) = 1 imply 

II Vporj ||o<|| Vporj ||3=|| Vp o ?7 o /3^^ ||3<|| Vpor] II3 (1+ || I3~^ \\^) < C. 
Proposition 4.16 gives 

II V^^^ o ^ llo < C II {D?^r^ II1+, (11 DrDf + II DlDf 

< C II ?? II2+. (II V/ II2+. + II V/ II 
Invoking (4.2) and corollary 4.19 we obtain the result. □ 

Finally we combine all of the above to obtain the desired energy estimate: 
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Proposition 4.22. 

Z II V/ \\l + II V/ 111 <C7 II V/ II? +^ f II V/ lloll V/ II2 
k 2 k Jq 



C 



+ II V/ lloll V/ ||2+e (II V/ 111 + II V/ h) (4.58) 

+ II V/ lloll V/ ||2+e (1 + k II V/ h+s +k II V/ ||3+e) 



Proof. From lemmas 4.13 and 4.12, propositions 4.17 and 4.21, and (1.12), it follows that 

1 II V/ \\l + II D^D.f \\l9< C II V/ II2+.II V/ 111 +C II V/ 11? +^ /* II V/ lloll V/ 
k 2 k Jq 



12 



+ J f W V/ lloll V/ ||2+e (II V/ 111 + 11 V/ 







Adding || V^f Hqq + || V/ 112+^11 V/ ||1 to both sides and invoking lemma 4.14 yields 

' 2 

\ II V/ \\l + II V/ ||1< C II V,/ WIq +C{1+ II V/ ||2+e) II V/ 1I2+.II V/ 111 

K 2 2 



+C II V/ llf +- II V/ lloll V/ II2 +- II V/ lloll V/ II2+. [\\ V/ 111 + II V/ II2 

+^ " "0" (1 + A: II V/ |l2+e +A: II V/ II3+.) 

The term || Vyf ||q q can be estimated by |1 ^uf \\1 g, which in turns is estimated by || V/ ||f because 
of (1.12). Since |1 V/ ||2+£ is small due to the bootstrap assumption, the term 

(1+ II V/ |l2+e) II V/ ||2+e|| V/ |l| 



can be absorbed into the term || V/ ||i on the left hand side, finishing the proof. □ 

2 

5. Regularity of V/ and elliptic estimates. 

In this section we shall prove propositions and 4.4 and 4.5. 

Proposition 5.1. Let r] and p be as in theorem 1.1. Suppose further that assumption (4-3a) holds 
for s < 41/2. Then V/ G H^{n). 

Proof. The strategy will be to solve for the highest derivatives of / in (4.25). Notice that /3 is as 
regular as rj. 

Recalling (4.12) and (4.13), 

D^Df = {Dlf + 2DrD,f - Drf) r + [DlD.f - 2Dlf - D,f) u, 

and so 

{DlDf, T + D^Df) = {Dlf + 2DrD,f - D^f) (l + D^f + D,f) 
+ {DlDJ - 2Dlf - D,f) {DrD.f - D^f) . 



(5.1) 
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Hence (4.25) gives 

{T,AoTjNo^ = -\D?j\-^ {Dlf + 2DrD,f - Drf) 

- \D?j\-^(^ [Dlf + 2DrD,f - Drf) (1 + Dlf + D,f) 
+ {DlD.f - 2Dlf - D,f) {DrDJ - Drf) 

- {DrD.f - Drf) ) (1 + Dlf + D,f) . 

We write this as 

(r, ^ o ^iV o 57) = -Dlf (\D^-^ + \D^-^ (1 + Dlf + D.ff 
- DlDJ (1 + Dlf + DJ) \D^~^ + 7^l, 
where IZi is a polynomial expression in at most two derivatives of /, and in \Dr]\~'^. Analogously, 
{v,AoljNo^ = {DlD.f - 2Dlf - DJ) + \D^-^ 

- \D?jr^ ( {Dlf + 2DrD,f - Drf) (1 + Dlf + D,f) 
+ {DlDJ - 2Dlf - D,f) {DrDJ - Drf) 

- {DrD.f - Drf) ) {DrDJ - Drf) , 

which can be written as 

{u, Ao^No?j) = -DlD.f (\D?j\-^ + \D?j\-^ {DrD.f - Drf)"" 

- Dl\D7i\~^ (1 + Dlf + DJ) {DrD.f - Drf) \D?i\-^ + 7^2 

where 7^2 is a polynomial expression in at most two derivatives of / and in \Drj\^'^. 

To get a formula for Norj we note that the unit tangent vector to drj{^l) is \Drr]\~^{T+DrV f). Thus 
if denotes a clockwise right-angle rotation of the plane we find thatiVo^= J(r + L»^V/) 

which equals IZ^^t/I-^z/ + JDrVf). Then letting W = JDrVf, we get 

N o?j = \DrV\~\i' + W). (5.3) 

Use (5.3) to combine (5.1) and (5.2) as 

ADlf + nDlDJ = -Ao?j{T,W)\Dr7j\'^ -TZi (5.4a) 
r2Dlf + BDlD^f = -Ao^DrVr^ - Ao^DrVr^ {i^,W)\DrV\'^ - 712, (5.4b) 

where ri and r2 are again polynomials in at most two derivatives of / and in \Drj\~^. Also A and B 
are given by: 

A = m-^ + (1 + Dlf + Dj)' , 

B = \D?}\-^ + \D?}\~^ {DrDJ - Drf)"" . 

Notice that W also contains at most two derivatives of /. 

Given assumption 4.1 for s < 41/2, if k is sufficiently large (and so V/ sufficiently small), use 
(4.27) to expand |-D??|"^ = 1 + Mi and |-Dr?|~^ = 1 + 2Mi + Mf, as it was done in proposition 4.15. 
It then follows that the lowest order terms in ri and r2 are 0{S/ f), and in particular ri and r2 can 
be made arbitrarily small by taking k large. It also follows that A and B are nowhere zero, so the 
algebraic system (5.4) can be solved for D^f and D^D^f . We conclude that 



(5.2) 
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r Dlf = ViiAo?j,{T,W),{i^,W),TZ3) (5.5a) 
I D^D.f = V2{Ao?},{T,W),{u,W),ni) (5.5b) 

where Vi and V2 are polynomials in their arguments and TZ^ and 7^4 are polynomial expressions 
containing at most two derivatives of / and in \DT-rj\~^. The explicit dependence on \DT-rj\~^ can be 
eliminated by expanding \DT-r]\^^ = 1 + M. Here M is a smooth function of D^f which is 0{D^f). 
Hence this term is absorbed into TZs and 7^4. 

We can now determine the regularity of /. First notice that /? is as regular as rj. By construction, 
we know that V/ is at least in H^^^^{Q) as well. From the hypotheses we have o ij £ H^{Q) and 
Vpo 07] e H^{i^), hence VAh o r? € H^{^), and V^j^ 07/ = VAh °i]o f3~^ G H^i^), where we recah 
that pq and Ah have been defined in (3.2). Since 

V{Ah orj) = VAh o rjDrj 

one sees that V{AH°'n) ^ so An^'f] G H^(n), and therefore Aorj G H^^/'^{dQ.). We infer from 

(5.5) — using also the fact that the right hand side of this expression has at most two derivatives of 
/ — that Drf G H^'/^{dn) and D^f G H^'/^{dn), which implies V/ E H^{n). □ 

proof of proposition 4-4- We know already that 77^ S H^^^^{il.) (see remark 1.2), so that V/^ and 
V/fc are in 

i/4V2(Q). As they satisfy V/fc(0) = = V/fc(0), shrinking Tk if necessary we can assume 
that assumption 4.1 holds for s < 41/2. We can therefore apply proposition 5.1 and conclude that in 
fact V/fc G H^{Q). Shrinking once more, we can now assume that assumption 4.1 holds in full, 
finishing the proof. □ 

proof of proposition 4-5: Either ^ = id +V/ or ^ = /3 + V/ o /3 = id + J^v o j3 + V/ o /3. || V/ Us will 
be small by the bootstrap assumption; for small time Jq || v o p \\^ will be small due to proposition 
4.17 and corollary 4.18, whereas || V/ o /3 II3 will be small in light of proposition 1.5, the bootstrap 
assumption and corollary 4.18. Therefore is near the identity in and so is a bounded elliptic 
operator from Hq{Q.) to H^~^{^}). for s < 3. □ 



Remark 5.2. Notice that in the previous proof is bounded uniformly in ^ as long as the bootstrap 
assumption holds. We shall show in the next section that equations (4.3) in fact hold in a uniform 
time interval, so the perturbation will be uniformly bounded in a small time interval independent of 
k. 



6. The bootstrap estimate. 

The goal of this section is to show that assumption 4.1 holds uniformly on a small time interval. 
As pointed out in remark 4.2, it is enough to our purposes to prove this for finitely many values of s. 

Proposition 6.1. Assume that (4.3a) holds for s = 0, 1, 2, 2 + e, 3 + e, 6, and (4-3b) holds for 
s = 0, 1, 41/2. Then there exist a number A > and a constant C > 0, independent of k, such that 

II V/fc ||,< -^—T, s = 0, 1, 2, 2 + e, 3 + e, 6, (6.1a) 

k 2 

II VA ||s< ^ = 0, 1, 41/2. (6.1b) 

k 2 ^'^ 
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proof for II V/ ||s, s = 0, 1, 2, 2 + e, 3 + e: We will use an iteration scheme, hence, write assumption 
(4.3a) as 

II V/ ||.< -^z^^ (6-2) 
A; 2 + 2 

for s = 0, 1, 2, 2 + e, 3 + e, 6, where is defined inductively by ^^^M = g and 



(5,+i(1.5) . r 0.5 + 5,(1) 1.5 + 5^(2 + e) -e 

— = mm < , 

2 I 2 ' 4 

5,+i(l) ^ 5,+i(1.5) 

2 3 
5,+i(2) _ 4<5,+i(1.5) 

2 9 
'^i+i(2 + £) ^ (4-£),5,+i(1.5) 

2 9 
5j+i(3 + £) _ (3-g)5,+i(1.5) 



}, (6.3a) 
(6.3b) 
(6.3c) 
(6.3d) 
(6.3e) 



Notice that with the exception of the first line, j ' + 1, and not j, appears on the right hand side of 
the above expressions. 

Plugging (6.2) into (4.58) yields 

rllv/||g + ||v/|||<-^+ ^ 



2 k^+SjW ,ia5+£i(2) 



k 2 



2 



C / 1 1 1 

+ 



lU.5-e I "jV^T'-; I , o I "J 1 e I ) 

2 ^ 2 y/v ^ i^^-t- 2 /j;2 2"'" 2 / 

For e sufficiently small, the term in parenthesis on the right hand side is bounded by a constant C. 
It then follows that the second term on the right hand side is bounded by the third term because 
from the above definitions one has + Sj(2+e) ^ lo^ _^ -^y^- Therefore 



1 11 „ A i|2 I II V7 ^ l|2 ' ^ ^ 



I II V/ llo + II V/ lis < -^^^ + ^ 

2 + 2 

c c 

+ 



(6.4) 



l.5+Sj{2+e)- 



;,4.5+0.5+5,(l) 5 
In particular, 

C C C 

II llf- 4.5 I + 4.5 I 1.5+5^{2+e)-e ^ 45 -Sj + rO^ " f^^-^) 

/C2"'"2 k 2 ^ 4 fc2''"2 

In order to handle the other values of s, recall the interpolation inequality 

II u \\\-''< C II u ff'W u II*-" for r<s<t, (6.6) 

which implies 

II V/ ||*-i-5< C II V/ lir^-^ll V/ 11*75^ for 1.5 < s < t < 6, (6.7) 

and 

II V/ ||l-^-''< C II V/ \\{~{\\ V/ for < r < s < 1.5. (6.8) 
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Using (6.2), (6.5), and (6.7) with t = 6 yields 

II V/ ii.< (6.9) 

A; 2 + 2 

for s = 2,2 + e, 3 + e, with 5(2), 6{2 + e) and (5(3 + e) given by (6.3c), (6.3d) and (6.3e), respectively. 
Then using (6.2) and (6.5) with (6.8) and letting r = yields 

II V/ ||i< -j^, (6.10) 
A;2+ 2 

with 5(1) given by (6.3b). 

Now one can repeat this argument: given a value of Sj{s) in (6.2), estimate (6.4) yields (6.5) with 
(5j-|_i(1.5) given by (6.3a); this then gives, through (6.9) and (6.10), newly improved values for the 
remaining 5j{s). One can check from the definitions (6.3) that it is possible to iterate this procedure 
a sufficient number of steps as to obtain (i) (6.1a) with some A', for s = 1, 1.5,2,2 + e,3 + e; (ii) 
^2+1^ > 1 + 7 and ^ > 1 + 7 for some 7 > 0. 

Now, since || V/ ||o<|| V/ ||i, (6.10) gives 

l|V/|lo< ^ 



giving the result with A = min{A', ^}. □ 

proof for II V/ ||s, s = 0, 1. Using estimate (6.4) and the above, 

Ck C 

II V/ ||o< 1 TTTT < — ■ (6.11) 

Using (6.11) with (6.8) and letting r = 0, s = 1 and t = 4.5 produces 

II lli^ , M + M^ ' 

giving the result with possibly a new constant A. □ 

Remark 6.2. With (6.1a) proven for s = 0, 1, 2, 2 + e, 3 + e and and (6.1b) proven s = 0, 1, one can 
apply the interpolation inequality (6.6) once more in order to obtain 

II V/fc |U< ^S-r, s < 6, (6.12a) 

k 2 

II VA ||s< -iSfTX^ ^ < 4-5 (6.12b) 

k 2 

possibly with new constants A and C. These can be assumed independent of s since we shall use 
(6.12) for only finitely many s-values. In particular these bounds are valid for integral s. We 
henceforth assume (6.12). 

proof for II V/ He- Prom (5.4a) 

^3^_ r^D^^D,f + Ao7j{T,W)\Dr7j\-'+n^ ^ ^^^^^^ 

From our constructions and application of (1.11) it is seen that for V/ sufficiently small 

II n hi/,,d + II (r, W) 1131/2,9 + II ^1 ||3V2,9< C II D^f ||3i/2,a< C II V/ ||4i/2,o< C II V/ lis (6.14) 



Drf \kl/2,d< C (II Drf h,d + II D'J ||3V,,s) < 7T^, 
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and 

II \DrV\-' hy2,d< C{1+ II D^f ||3V2,a) < C(l+ II V/ lis) < C. (6.15) 
Furthermore, from proposition 4.15 with s = 3 one obtains 

M o ^ \y/,^Q< C(l+ II V/ ||5v,,a) < C(l+ II V/ lie) < C, (6.16) 
provided that V/ is sufficiently small. We also have 

II D^rDuf ||3i/2,a< C II V/ ||5V2,a< C || V/ ||6< C. (6.17) 
Using (6.14), (6.15), (6.16) and (6.17) along with (6.12) and (6.13) yields 

II D'rf ll3V2,a< (6-18) 

for some A > 0. 

We shall now use equation (2.2a) to estimate /. Since this equation is unchanged by subtracting 
a constant from /, and since our goal is to estimate V/, we can assume that Jq^ / = 0. Then 
Poincare's inequality implies 

II / 1161/2,9- II Drf ||5i/2,9 • 

But 

C 
so that 

II / ll6V2,9< (6-19) 
where (6.12) and (6.18) have been used. From (2.2a), 

II V/ lie < C II / ||7< C( II A/ lis + II / ||6v,,a 

<C[\\ {D^ff 115 + 11/1161/2,9 

In order to estimate || (D'^f)'^ Us, let us look at the highest order term 

II D'^D'ff llo < II {D'D^f){D'~'D'f) \\, 

1=0 
4 

<CY,\\ {D'+'Df){D'-'Df) llo 
e=o 

Applying (1.11) with s = l+e, r = and the term with the r-norm on the the right of (1.11) being 
always the one with fewer derivatives in the product {D^~^^ D f){D^^^ D f) , leads to 

2 

II D'iD^ff ||o<C^|| V/ II 2+^+ell V/ lle-^ 

<\\ V/ ||3+e|| V/ lis + II V/ ||4+e|| V/ ||4 + || V/ ||2+.|| V/ He 

^7^+ II W l|2+e|| V/ lie. 



(6.20) 
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The remaining terms in || (D'^f)'^ \\^ are estimated similarly, yielding 



{D^ff h< -S-+ II V/ II2+.II V/ lie . (6.21) 



Combining (6.19) (6.20) (6.21) gives 



II V/ ||6< II V/ II2+.II V/ lie 

The term || V/ ||2+e|| V/ He can be absorbed on the left hand side by (6.12a), finishing the proof. □ 
proof for II V/ ||4i/2. Differentiate (6.13) in time and proceed as in the proof of || V/ \\q. 

7. Proof of theorem 1.1. 

We are now in a position to prove theorem 1.1, assume therefore all its hypotheses. As we are in- 
terested in passing to the limit, it will be convenient to reinstate the subscript k ^. The t-dependence 
will also be written when we want to stress the interval where certain estimate is valid. We shall 
make extensive use of proposition 1.5, so we remind the reader that J{Pk) = 1- 

proof of theorem 1.1: For each /c, let Tk be the maximum interval of existence for (1.1). Since 
((0) = id = r/fc(0) and 77(0) = m = rik, 

II C{t) - Vk{t) II41/2 =11 m - C(0) +?/fc(0) - r/fc(t) II41/2 

<ll C{t) - C(0) II41/2 + II Vk{t) - r/fc(0) II41/2, 

and 

II at) - 77fc(t) 1141/, =11 c(t) - c(o) + %(o) - f]k{t) 1141/2 

<ll C(t) - C(0) II41/2 + II %(i) - %(0) II41/2 
As C; C) Vk and fjk are continuous functions of t, given > we can choose = Tkig) such that 

II C{t)-Vk{t) ||4i/2< Q, on [0,Tfc], (7.1) 

and 

II C{t)-Vkit) ||4i/2< Q, on [0,Tfc]. (7.2) 
Since C exists for all time, there exists a constant M such that 

II C{t) II41/2 + II Ut) ll4i/2< M, for t < Tk (7.3) 
Combining (7.1), (7.2) and (7.3) yields 

II m{t) II41/2 + II m{t) ||4i/2< 2M + 2^ = M„ on [0,Tfc]. (7.4) 

By proposition 4.4, the equations (4.3) hold on a time interval [0,7fc], and shrinking Tk if necessary 
we can assume Tk < Tk- Proposition 6.1 then implies that equations (6.1) hold on Tk- Therefore 
from (3.1) , (7.4) and proposition 1.5 it follows 

II /3fc(t) ||4i/2< M, + C\\ Vfk{t) II41/2 (1+ II fik{t) II41/2), on [{),Tk]. 



Notice that the operators Li and L2 also depend on fc, but we shall not attach a subscript k to them as to avoid 
a cumbersome notation. All we shall need about Li andL2, in fact, is that they are bounded operators with norm 
bounded above by a constant independent of fc — a condition that clearly holds for k large. 
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Because of (6.1a), if k is sufficiently large, the term C \\ V fkif) II41/2 || /3fc(f) II41/2 can be absorbed on 
the left hand side and therefore 

II Pk{t) \W/2< CM, + C II Vfk{t) ||4i/2< C, on [0,rfe], (7.5) 

where (6.1a) has been used to bound || Vfk{t) ||4i/2- Then, using (7.5) along with proposition 1.5 
and (4.2), (3.5), (4.19), (1.11), (6.1b) gives 

II Vk{t) II41/2 <|| rjk{t) o Pk^{t) II41/2 + II D^^(t)Vfk{t) II41/2 + II V/fc(t) II41/2 
< C+ II Vk{t) II41/2 II Vfk{t) II51/2, on [CTfe]. 

As before, the last term can be absorbed on the left hand side if k is large enough because of (6.1a), 
so that 

II Vk{t) II41/2 < C, on [0,71^]. (7.6) 
From (7.5), (7.6) and proposition 1.5 it immediately follows that 

II /3fc(t) ||4i/2<C, onlCTfc]. (7.7) 
Therefore, invoking again (3.1) we obtain, with the help of proposition 1.5, (7.5) and (6.1a), 

II ilk{t) - m) ll4V2< C II Vfk{t) ||4i/2< ^, on [0,rfc], (7.8) 

and from (3.5), (7.6), (1.11), proposition 1.5, (6.1a) and (6.1b) we have, after recalling that I3k{t) = 
Vk{t)opk{t), 

II %(t) - h{t) II41/2 < C II V/fc(t) II51/2 +C II Vfk{t) ||4i/2< ^, on [CTfc], (7.9) 

where A comes from proposition 6.1. From (7.8) and (7.9) we see that in order to obtain the desired 
convergence it will suffice to show that Pk{t) C{t) a-iid /3fc(i) C,{t)- We proceed therefore to 
analyze the equation of motion for /3k{t), i.e., (3.9b). Differentiating /3fc = f o in time, 

/3k = {vk + yv,Vk)o(3k. (7.10) 

But since 

Vk + P{^v,Vk) = Vk + {I - Q){^v,Vk) 

composing (3.9b) with (3^ on the right and using (7.10) yields 

f3k = (Q(V^,,„^-i A o /3-^)) o /3fc - 2{L^'PV,,Vfk) o (3k - [L^^PD\„^V fk) o (7.11) 

We recall that the pressure p appearing in (1.2a) is determined by the fluid velocity § = C,o C~^- 
fact, exactly as in (7.10), we have C = (^^ + V^t?) o so that (1.2a) becomes 

,? + Vtfi? = -Vp. (7.12) 

Taking divergence of (7.12), using div(t?) = and recalling the definition of Q, (1.2a) can be written 
as 

c = (Q(v^oc-iCor'))oc. (7.13) 

Furthermore, since Pk ° Pk^ and ( o are on the image of P, we see that (7.11) and (7.13) can also 
be written as 

(3k = Z{PkJk) + nk, (7.14) 
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and 

C = Z{C,0, (7.15) 

where Z is the operator defined by 

Z{a,v) = (Q{S/^oa-^P{v o a"^))^ °a. 

and TZk is given by 

7^fc(^) = -2(L-1PV„,VA) o/3fc - {L-'PDl,^Vfk)o(3f,. 

Z is well defined since if w is a H'^ divergence free vector field, then the gradient part of ^wW, i.e., 
QCVww) is also in and not merely in H^~^ as one would expect at first sight. This follows from 
the definition of Q since the divergence free condition implies that div(V^^i') depends on only first 
derivatives of w; see [22] for details. 

Integrating equations (7.14) and (7.15), subtracting, and using that /3a:(0) = uq = C(0) (because 
V/fc(0) = 0) gives 

II $k{t)-C{t) 1141/2 < / II Z{^kJk){s)-Z{C,0{s) II41/2 ds+ [ II TZkis) 1141/2 ds, (7.16) 
Jo Jo 

We need to show that |[ TZ^ II41/2 is small. For this, first notice that 

PD.,Vh = [P, A-JVA + A-,PV/ = [P, A-JVA, (7.17) 

where [■, •] is the commutator and we used that P o V = 0. Recall that if £1 and C2 are differential 
operators of orders ii and £2, respectively, then their commutator [i2i,>C2] has order ii + £2 — 1- 
Therefore [P, A*.] is a zeroth order operator and hence (7.17) can be estimated in if^V2 with the 
help of (1.11): 

II PA.VA ||4V2< C II Vk II41/2II V/fc II41/2 . (7.18) 
Analogously, [P, Afew^] is first order and therefore 

II PP'^kVk II41/2— II Vk II41/2II V/fc II51/2, (7.19) 

where we used again (1.11) and the fact that the operator A^t^fc does not depend on derivatives of 
Vk (see (3.8)). Combining (7.18) and (7.19) with (7.5), (7.6), (6.1), and proposition 1.5 we obtain 

II 7^fc(^) ||4i/2< ^, on [CTfc]. (7.20) 
Then, from (7.20) and (7.16), it follows 

II 4(t) - m II41/2 + I II ^(A> A)(s) - Z{C,C){s) II41/2 ds, on [0,Tk]. (7.21) 

In [16] is proven that the map Z is smooth. Since 

(C(i) , C(i) ) ( [0, Tfc]) is compact in P^'^' {ft) x H^'/^ (n) , 
the map Z is uniformly Lipschitz in a small neighborhood of the curve (^(t), C(i)). Therefore, because 
of (7.1), (7.2), (7.8) and (7.9), if g and Tk are sufficiently small and k sufficiently large, 

II ziPkJkM - z{cxm 1141/2 < c II (3k{t) - m 1141/2 +c II m) - m lUv^ on [o,7X]. (7.22) 

But, 

II Mt) - m II41/2 < / II $k{s) - C{s) II41/2 ds. (7.23) 
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Using (7.22) and (7.23) into (7.21) we obtain, after iterating the inequality (7.21), that 

II Pk{t) - m + II $k{t) - at) lUi/, < ^e^\ on [0, Tk]. (7.24) 

Recall now that the time interval was defined by the validity of (4.3). Proposition 6.1 shows that 
this implies that a better estimate, namely, (6.1), holds on [0,7fc], and thus if k is sufficiently large, 
7fe is a non-decreasing function of k, as long the solution %(t) exists on that interval. But now, since 
the solution exists for all time, it follows from (7.24), equations (7.15), (7.14) and (7.20), that the 
interval of existence of /?fc(t) is also non-decreasing with k (for k large), and then the same is true for 
Vkit) by (7.8) and (7.9). We conclude that there exists a T > 0, independent of k, such that r]k{t) 
exists on [0, T] for all k, and such that (after combining (7.24) with (7.8) and (7.9)), 

sup II r]k{t) - C(t) II41/2 + sup II f]kit) - C{t) II41/2 — > as k ^ 00, 
0<t<T 0<t<T 

finishing the proof. □ 

7.1. Final remarks. Here we expand on several comments made in the introduction, and make 
some further remarks. 

We start by pointing out how the theorem can be generalized to include different initial conditions 
UQk, not necessarily tangent to the boundary, but which converge suitably to an initial velocity uq of 
problem (1.2). In this case, we would have 

Uok = P{uok) + Q{uok) = Vok + V/ofc. 

We then have to choose 

vok - ""0 -p- and Vfok ~ -p-, 

in the appropriate norms. It is not difficult to see then that if A is large, the estimates we presented 
will still go through. In fact, in the energy estimates, there will be an additional term uq^ coming from 
V/fc(0), which can then be grouped with the other terms in ^ without affecting the result, provided 
that A is suitably chosen for this purpose. Similarly, when we considered || ({t) — r]k(t) II41/2 and 
|[ C(*)~/3fc(i) 1 1 41/2 ill the previous proof, we would have extra terms of the form || UQ — {vok+^fok) 1 1 41/2 
and II tio — 1 1 41/21 respectively which again will be negligible errors of the order 

The importance of having only one boundary component with constant mean curvature at time 
zero can be seen from the boundary condition for the pressure p. p is explicitly coupled to rj via 
its gradient. The "boundary part" of p is given by the harmonic extension Ah of the the mean 
curvature A, and therefore "VAh = when A is constant so one can expect VAh to be small for 
small time. 

We point out that our method can also be used to prove well-posedness in £"^{0,). This is done by 
using an iteration scheme for the system (3.9). The iteration will converge if / is sufficiently small, 
which, due to the estimates of section 2, will be the case provided that Sq in theorem 2.1 is small 
enough. This will be shown in a future work [21]. 

We finally comment on the higher dimensional case. There seems to be no conceptual problem in 
generalizing our proofs to three and higher dimensions, the difficulties being mostly computational. 
Notice that in section 7 we did use the fact that solutions to (1.2) exist for all time, a feature unique 
to those equations in two spatial dimensions, but this is superfluous: all we need there is a fixed 
time interval [0, T,,] in which (1.2) has a unique solution (see [16, 19, 20] where a similar convergence 
argument is used). In fact, we expect that most of the arguments here presented can be generalized 
to higher dimensions in a more or less straightforward manner. The part that presents substantial 
difficulties is the derivation of the energy estimates of section 4, in particular those dealing with the 
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mean curvature of the boundary. The simphcity of the mean curvature operator for one-dimensional 
boundaries was employed to a great extent in section 4. The expressions will become much harder 
to handle in three and higher dimensions, as it should be expected from usual difficulties posed by 
mean curvature type of equations. 
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